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ABSTRACT 

We present an analysis of our Chandra Low Energy Transmission Grating 
Spectrometer (LETGS) observation of the quasar MR 2251 — 178. The warm 
absorber of MR 2251—178 is well described by a hydrogen column density, 
]Vh ~ 2 x 10 21 cm -2 , and an ionization parameter log(£) ~ 0.6. We find in 
the spectrum weak evidence for narrow absorption lines from Carbon and Nitro- 
gen which indicate that the ionized material is in outflow. We note changes (in 
time) of the absorption structure in the band (0.6 — 1) keV (around the UTAs 
plus the O VII and O Vlll K-edges) at different periods of the observation. We 
measure a (0.1 —2) keV flux of 2.58 xl0~ n ergs cm~ 2 s _1 . This flux implies that 
the nuclear source of MR 2251—178 is in a relatively low state. No significant 
variability is seen in the light curve. We do not find evidence for an extra cold 
material in the line of sight, and set an upper limit of Ah ~ 1.2 x 10 20 cm -2 . The 
X-ray spectrum does not appear to show evidence for dusty material, though an 
upper limit in the neutral carbon and oxygen column densities can only be set 
to N c l rj 2 x 10 19 cm -2 and N Q l « 9 x 10 19 cm -2 , respectively. 

Subject headings: galaxies: active - X-rays: galaxies - quasars: individual 
(MR 2251-178) 



1. Introduction 



Warm absorb ers have provid ed deep insights into the nuclear environment of Active 
Galaxies (AGNs). iHalpernl (119841 ) reported its presence for the first time using the Einstein 
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observation of the QSO MR 2251 — 178. Since then, warm absorbers have been commonly 
found in about 50 % of the AGN spectra (see iKomossa fc Hasingerl 120031 ; ICrenshaw et al. 



20031 . for reviews), and their study has enriched our knowledge about the ionization and 
the kinematics of the gas composing these systems, important for the understanding of the 
evolution of these objects, and the AGN unification picture. 



In this context MR 2251-178, at z = 0.06398 (Bergeron et al.lll983 ). has been the 
subject of various stud ies, first due to its historical importance as the first qu asar discovere d 
by X-ray observations (IRicker et al.lll978l ). the first warm absorber reported (lHalpernlll984l ). 
and because it displays a number of outstanding characteristics. It is surrou nded by a large 



[O il l] emission line region, located in the outskirts of a cluster of galaxies (IBergeron et al. 



19831 ). with a high ratio L x /L opt . Sev eral studies have been focused on the characterization 
of the absorbing material properties. iKomossal (120011 ) fitted a warm absorber model to the 
ROSAT observation of MR 2251 — 178, and found a high ionization state of the absorber 



(log U ~ 0.5), and a column density of Ah ~ 10 22 6 cm 2 . iMonier et al.l (120011 ) reported the 



pres ence of absorption line s of N v A 1240 and C IV A 1549 blueshifted by few hundreds of km 



Ganguly et al.l (120011 ) noticed variability of the C IV A1549 line by taking observations 



from the Hubble Space Telescope (HST) 4 years apart, and inferred a maximum distance 
from the source of ~ 2.4 kpc. 



More recently, iKaspi et al.l (120041 ) established a scenario where clouds crossing our line 
of sight show up by the presence of absorption lines coming from high ionization states, 
displaying a wide velocity range from to ~ 600 km s -1 , taking a 8.5 yr set of data from 
ASCA, BeppoSAX and XMM-Newton. 



Gibson et al.l (120051 ) analyzed data from the High Energy Transmission Grating on board 



Chandra, and found evidence in the spectrum of MR 2251—178 of highly ionized, high- 
velocity (v ~ 12, 00 — 17, 000 km s" 1 ) outflowing material. They report similar column 



densities to those of IKaspi et al.l (120041 ). of a few 10 21 cm 2 , and establish conditions for the 



accretion and mass-loss rates. 



From the Fe xxvi L Q line, iGibson et al.l (120051 ) concluded that, unless the absorber cov- 
ering factor is very low, the mass-loss rate is approximately an order of magnitude higher than 
the accretion rate that would provide the radiation power of MR 2251 —178 (~ 0.2 M G) yr ~ ] 



at 10 % efficiency). This is a different conclusion from that reached by lMonier et al.l (120011 ). 
in the Hubble Space Telescope (HST) analysis of the same object, where they established that 
the accretion and mass-loss rates are essentially the same. The inquiry into this point is of 
relevance to the understanding of the evolution mechanisms in the nucleus of MR 2251 — 178. 
On the other hand, the HETGS high resolution observation also reveals the presence of 
emission lines coming from highly-ionized species, and shows that the emitting material is 
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not in our line of sight, opening the possibility that the absorption and the emission lines 
have different physical origins or that they are located at different distances from the nu- 
cleus. This phenomeno n has been observed in NGC 4051 in the optical/ultraviolet band. 
Komossa fe Finkl (119971 ) concluded that the coronal emission lines from lowly ionized species 
of Fe ([Fe vn]-[Fe xi]) observed in the spectrum of this object, have a different physical ori- 
gin from the warm absorber, consistent with the fact that they could come from a different 
spatial regio n, the warm absorb er located at the nucleus and the coronal lines extended out 
of ~ 150 pc JNaeao et al.lboooh . 



Finally dusty warm absorbers have been found in a number of AGN (IKomossa fe Hasinger 



20031 . and references therein for a review) and are predicted to l eave their mark on the 
X-ray spectra of AGN , for instances the K- edges of O I an d C I ( IKomossa fe Finkl 119971 ; 
Komossa fe Baddll998l ). and the Fe-L edge (ILee et al.ll200ll ). For the first time, we search 
for the presence of a dusty warm absorber in MR 2251 — 178. Based on ROSAT data, IKomossa 
( ]200ll ) speculated about the presence of a second dusty warm absorber in MR 2251 — 178, 
this second warm absorber dusty, but the data did not allow multi-component fitting, which 
would provide important implications on the global structure of MR 2251—178. 

It is clear from this introduction that each of the previous missions and the use of more 
powerful spectroscopic instruments lead to new perspectives regarding the properties and 
the evolution of the nuclear environment of MR 2251 — 178. We present an overview and an 
analysis of the spectrum of MR 2251 — 178 as it is seen with the Low Energy Transmission 
Grating Spectrometer (LETGS) on board Chandra. The spectral resolution power of the 
instrument allows us to confirm some of the previous conclusions about the ionization and 
the kinematics of the warm absorber of MR 2251 — 178, and shed new light on the evolution, 
ionization and composition of this system. 

Technical details on the observation and data reduction are given in §2J Our spectral 
analysis is given in §3.21 The mean properties of the warm absorber of MR 2251 — 178 are 
presented in §3.4| with the corresponding spectral absorption lines analysis. To finish, we 
made a temporal analysis on the absorption structure around the O VII and O vm K-edges 
in §U We discuss the results in §5] and conclude in §HJ Throughout this paper, we use the 
following cosmological parameters: H = 70 km s -1 Mpc" 1 , Qm — 0.3 and Q\ = 0.7. 



2. Observation and Data Reduction 



We obtained a ~ 80 ks exposure time observation of M R 2251 — 178, performe d with the 
Low Energy Transmission Grating Spectrometer (LETGS, iBrinkman et al.ll2000l ) on board 
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Chandra (under the sequence number 700405 and OBS ID 02966). A log of the observation 
is presented in Table [TJ The spectrum was obtained by reducing the data with the Chandra 
Interactive Analysis of Observation (CIAC0) version 6.6. 

We use the default spectral extraction region (i.e., a "bow-tie" shaped region). When the 
LETG {Low Energy Transmission Grating) is used with the HRC-S detector [High Resolution 
Camera), this comprises a central rectangle abutted to outer regions whose widths increase as 
the dispersion distance increases. The background region is taken from above and below the 
dispersed spectra. The region shape for both, source and background negative and positive 
orders is precisely given in the file letgD1999-07-22regN0002 . f its, that can be found at 
the LETGS calibration webpage B 

Having properly extracted both, source and background spectra from each arm, we 
merged them, obtaining added source and background spectra. This is intended to increase 
the signal-to-noise ratio (S/N) of the final spectrum, and the spectral analysis (throughout 
this work) is based on this co-added spectrum. The effective areas (EA) for orders from 2 to 
10 used in fitting procedures were taken from the LETG+HRC-S effec tive areas webpage 2 . 



For the first order we used the corrected EA of iBeuermann et al.l (120061 ). EA orders from 1 to 



10 (positive and negative orders) were summed to be used with the corresponding co-added 
spectrum. 

The nominal LETGS wavelength range is 1.24 — 175 A (0.07 — 10 keV). However, we 
restrict our spectral analysis to the range 1.24 - 124 A (~ 0.1 - 10 keV). In the range ~ 
0.07—0.09 keV the S/N ratio is low (« 1.4 at ~ 0.09 keV), and we exclude bins in this band 
from our analysis. Also, we exclude bins that fall in the bands where the gaps of the detector 
(HRC-S) are located (i.e., 52—56 A and 62—66 A, for left and right gaps respectively). Figure 
[U shows the summed (positive and negative dispersion orders) background subtracted count 
rate spectrum of MR 2251 — 178 adaptively binned to have at least 100 counts from the source 
per bin over the range considered in this work. It is known that an instrumental feature 
~ 2.08 keV, related to the mirror, could be modeled by an edge with negative optical depth 
%. We model this Ir M edge of the mirror by using an edge model with energy fixed to 2.08 
keV, and an optical depth equal to —0.15. So all our fits include this negative edge model. 
All wavelengths and energies are presented in the observed frame. 



^ttp : //cxc . harvard . edu/ciao 
^http : //cxc . harvard . edu/cal/Letg 



http : //cxc . harvard . edu/ ccw/proceedings/03_proc/presentations/marshall2/s001 . html 



- 5 - 



2.1. Source Extent and Light Curve 

In order to search for evidence of a spatial extent of the source or the presence of other 
possible X-ray sources, for example jets, close to MR 2251—178, we take the zero-order image 
(Chip 2 on the HRC-S) of MR 2251-178 centered at RA= 22:54:05.08 and Dec(J2000)= - 
17°34'55" in a field of view of 10 "x 10 ". The extent of the emission (and the 90 % encircled 
energy; computed with the celldetect command of CIAO0) is consistent with being inside 
the point-spread function (PSF) dominated zone (for the LETGS ~ 1.8 ). We conclude that 
no evidence of appreciable extent of the source is found in this observation (we also checked 
for angular dependence in the X-ray image, and found that the emission of the object is 
symmetrically distributed) . 

The light curve constructed from the zero-order of the observation is shown in Figure 
El (top panel). Counts are binned to 200 s bins. We find a mean value of 111 counts bin -1 
(horizontal red line) without any significant variation within la (horizontal green lines). In 
the bottom panel we present the light curve of the background (see definition below), which 
shows that strong count variations occur in several time intervals. Data from these intervals 
were excluded from the spectral analysis presented below. We removed data coming from 
time intervals whose count rate were above 3cr (horizontal green lines) from the mean value 
(horizontal red line). The final spectrum was constructed from these good time intervals 
(GTI), time intervals where counts are within 3cr from the mean value of the fluctuation 
exhibited by the background. A circular region around the zero-order image with a radius of 
10" was used to extract the source light curve. The background extraction region is defined 
by four circles with radii of 25" equidistantly distributed from the source in order to avoid 
contamination from events associated with dispersion and cross-dispersion directions. 



3. Global Spectral Models 

3.1. Flux and Luminosity (0.1 — 10 keV) 

To compute the hard, and soft X-ray fluxes we fit a simple power-law to the total spec- 
trun@ absorbed only by the Galactic column density of A^h = 2.77xl0 20 cm -2 ( Lockman &: Savage 



http : / / cxc . harvard . edu/ ciao3 . 3/ download/ doc/ detect_mamial/cell_run . html 



5 The added spectrum described in section [5] adaptively binned to have at least 100 counts from the source 
per bin. The fit was carried out in XSPEC 11, and the fluxes and luminosities computed with the commands 
flux and lumin respectively (the used band was 0.1 — 10 keV, and the model includes the Ir edge at 2.08 
keV with r = -0.15). We note that the LETGS spectrum of MR 2251-178 is dominated by the soft X-rays 
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19951 ). and we obtain a photon index of T x = 1.94lo'o3 an d a normalization value of 5.64 



Iqos x 10 3 photons keV 1 cm 2 s 1 at 1 keV. Then we integrate the unabsorbed power- law 
| to obtain a flux in the 0.1 - 2 keV band / (ai _ 2 kc v) = 2 -58±g;g| x 10~ n ergs cm" 2 s -1 , 
and the (0.1 — 2) keV luminosity is L(o.i-2 keV) = 2.54 x 10 44 erg s -1 , a factor of ~ 1.5 
fainter than the ROSAT soft X-ray observation. The (2 — 10) keV flux is /(2-10 keV) — 1-64 
-n ils x ]_q-ii ergs cm -2 s -i ) anc [ a (2 — 10) keV source luminosity L(2-io kov) = 1.62 x 10 44 



_Q 05 a iu eigs cm s , <±iiu ci — iuj kcv suuice luminosity 1^(2-10 koV) 
erg s~ l . This is ~ 30 % lower than the 2 — 10 keV luminosity measured by iGibson et al. 
( 120051 ) of iv( 2 _io keV) — 2.41 x 10 44 erg s _1 (but see footnote 5). In our computation, 284 of 
the total of 321 bins included in this analysis are in the 0.1 — 2 keV band, so all the spectral 
analysis is actually driven by the soft part of the spectrum. The information contained in 
the 0.1 — 0.5 keV band is important to conclude that the soft X-ray flux is between ~ 10 — 80 
% of the hard fluxes historically measured from the spectrum of MR 2251—178. Recently, 
this object has been observed in the hard X-ray band (14 —195 keV) with the Swi ft Burst 
Alert Telescope (BAT), with a L(i4_igs k c v) — 10 45 erg s" 1 (IMarkwardt et all 120051 ). Taking 
our best-fit power-law and integrating over the 14—195 keV band, we predict a luminosity 
£(14-195 keV) = 3 x 10 44 erg s _1 , a few factor lower than the Swift observation. The full 
range (0.1 — 10 keV) luminosity, is £(0.1-10 keV) = 4.17 x 10 44 erg s -1 . Our measured (0.1 — 2) 
keV flux allows us to conclude that MR 2251—178 is in a relatively low soft X-ray state. 



3.2. Spectral Models 



The rest of the spectral analysis of the present work was carried out with the Interactive 
Spectral Interpretation System package (ISIS ver. 1.4.4) using the XSPEC modules from 
the lheasoft libraries^. 

MR 2251—178 is surrounded by a giant emis si on-line nebula o f ~ 10 — 200 kpc extent 



|e.g., iBergeron et al.lll983t iMacchetto et al.lll990l ). iBergeron et al.l (119831 ) estimated a col- 
umn density of a few times 10 21 cm -2 of the emission-line gas. They speculated about the 
presence of an extended neutral H I halo around MR 2251 — 178 which should be detectable 
in H I 21cm observations, or in the soft X-ray band, if also located along our line-of-sight 



band (i.c, < 2 keV). So there might be an underestimate of the flux in the band 2—10 keV. The error can 
be estimated around 20 %. 

6 We also compute the flux using a two powerlaws model, which represents in a better way the 2 — 10 keV 
band, with deviations from the data not larger than 10%. The resulted fluxes do not change significatively. 

7 http : //space . mit . edu/CXC/ISIS/index . html 

6 http : //heasarc . nasa . gov/ docs/sof tware/lheasof t/ 
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(see also Sect 6.1 of iKomossal 120011 . for further motivation). We begin our fitting procedure 
by searching for, and placing constraints on, any cold material along our line-of-sight to the 
quasar which is in excess to the Galactic absorption toward MR 2251—178. 

We fit a simple power-law, absorbed by a column of cold gas in the line of sight, to 
the X-ray spectrum of MR 2251 — 178. Both, the power-law parameters and the column 
density are free parameters. The best-fit photon index and column density we obtain are 
r x = l-79±g;g| and A" H = 1.68±g;g? x 10 20 cm" 2 , respectively. We plot in Figure El the best-fit 
absorbed power-law model along with the spectrum of MR 2251 — 178 in energy space. The 
best-fit column de nsity is about 40 % lower than the Galactic value toward MR 2251—178 
(2.77 xlO 20 cm -2 , lLockman fc Savage! Il995l ). If we fix the column density to the Galactic 



value and re- fit, we obtain a slight change in the value of the photon index to T x ~ 1.9 and 
the fit gets worse w ith a change in A% 2 ~ 140. Taking the upper limit of the extra cold 



material reported by IKomossal (120011 ) of AA^h = 5 x 10 19 cm -2 and including it in the column 
of gas, we obtain a photon index slightly steeper and the fit gets worse if we increase the 
column density of the extra cold material by one and two times AiVn (see rows 3 and 4 in 
Table [2]), increasing y 2 (respect with A^r- free to vary) by A% 2 = 260 and 396 respectively. 

We also carefully investigate if excess X-ray absorption, above the Galactic value, and a 
blackbody component could be combined to compensate each other so as to mimic a less ab- 
sorbed spectrum. This case is shown graphically in Figure [3] in solid for the composite model 
(i.e., zphabs x [pl+bb] in comparison with the absorbed power-law (top panel). The 
best-fit intrinsic (AGN frame) extra column density has a value (given the data) consistent 
with zero. 

We tried a last model composed of two power-laws and cold absorption. The result 
is shown in Figure [3] (second panel from top to bottom). A steep power-law (T x ps 3.5), 
with help of the cold absorption, describes relatively well the soft band of the spectrum 
(~ 0.1 — 0.6 keV). A flatter second power-law with T x ~ 1 is responsible for describing the 
hard band (~ 0.7— 10 keV). An extra absorption, above the Galactic value, arises naturally 
from the fit (with the same number of parameters as the pl+bb model). A material with 
column density A^h ~ 1.2 x 10 20 cm -2 accounts for this extra absorption (see the best-fit 
column density in Table [2] row 6). We enforced this model to have an excess absorption of 
this best-fit column density plus 1 and 2 times A Ah, and re-fit, obtaining worse fits (see rows 
7 and 8 in Table [2]). We conclude that there is little evidence for cold excess absorption with 
a column density of A^h ~ 10 21 cm -2 . The maximum amount of cold absorption we can hide 



9 power-law plus blackbody modified by a column of gas in the AGN-frame. The model is corrected by 
the Galactic absorption. 
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in the X-ray spectrum of MR 2251 — 178, based on the model fit involving two powerlaws, is 
iVn ~ 1.2 x 10 20 cm -2 , and we consider this as a safe upper limit to the excess absorption 
along the line of sight at the epoch of observation. 

From now on, all the models considered in this study are multiplied by an absorption 
column of gas (phabs model in XSPEC0) with = 2.77 x 10 20 cm -2 to describe the 



Galactic absorption toward MR 2251—178 (ILockman fc Savage! 1 1 9951 ). First, we focus our 



attention on the fit of a power-law to describe the 1.24—124 A band of the spectrum. A 
single power-law gives T x = 1.94+^2, with large residuals (i.e., 15—40 %) in the ~ 1.2—7 A 
and close to 18 A (see Figure HJ), and a poor statistical fit quality of xl — 3.00. We included 
a thermal component to account for any soft X-ray excess. A blackbody with a temperature 
of kT = 81.5^2 e V significantly improved the fit with xl — 1-63 (for d.o.f = 305). With 
this addition, the slope of the power-law changes to a flatter value of T x = 1.53lo!o2' anc ^ the 
model agrees with the data within 30 % almost from ~ 1.2 to 40 A, except for the residual 
around 18 A. 



Now, let us consider the absorption features at » 18 A and ~ 15 A in turn. iKomossa 



( 120011 ) fitted a power-law with two edges at the theoretical positions of O VII and O Vlll, 
739 eV and 871 eV respectively, representing the warm absorber, obtaining Tq V n — 0.22 ± 
0.11 and r viii = 0.24 ± 0.12. The high sensitivity and resolution of the Chandra LETG 
spectrometer coupled with the good quality of the observation of MR 2251 — 178, allow us 
to fit these features with a high precision. Our plx 2 edges (« 18 A and ~ 15 A in the 
observed frame) model gives us a significant improvement (xl = 1-38 for 305 d.o.f) over the 
fit of the power-law alone (xl — 3.00 for 307 d.o.f) and the power-law plus the blackbody 
component as well (xl = 1-63 for 305 d.o.f). This model is in agreement with the data ~ ± 
15 % (data/model ~ 1.00 ± 0.15, Figure HJ) over almost the entire spectrum, except for the 
wavelength range 1.2—6 A (and possibly at wavelengths > 60 A , but see the large errors). 
The continuum optical depths are r Q VII = 0.56+q;q4 and r Q vm = 0.22+0^. We compare these 
values with those obtained from ROSAT. In the case of O Vlll they appear to agree within 
the error bars. However, in the case of O VII the optical depth is approximately a factor of 
2.5 higher. In this direct comparison model-to- model, the Chandra data suggests that the 
contribution of this feature has changed over the years (~ fifteen). 

Finally, we can build a model combining all the ingredients mentioned above leading us 
to a global view of the spectrum: A power-law plus a soft-thermal component modified by 
two warm absorber edges. This model is shown in Figure H] and the maximum deviation of 



lc http : / /he asarc . nasa . gov/ docs/xanadu/ xspec/manual- The abundances are solar from 
Anders fc Grevessel |l9891 ) 
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the ratio data/model is ± 15 % over almost the entire spectrum, with xl — 1-25 (for 303 
d.o.f). The F-test identifies this improvement (the inclusion of the black body emission) 
as significant (at > 99% significance level), over the plx 2 edges model. The blackbody 
has a temperature of kT ~ 50 eV with the slope of the power-law being in agreement with 
other models of this study (~ 1.7 — 2) and with slopes reported by p revious missions (i.e . 
ROSAT), but not with the flatter slope of ~ 1.4 reported recently by iGibson et al.l (120051 ) 
in the HETGS observation of the same object. We can see that the inclusion of the strong 
O VII feature displaces the temperature of the blackbody from kT ~ 80 to 50 eV with respect 
to the pl+bb without edges. The reduced temperature of the bb helps to improve the fit in 
the range ~ 20 to 60 A, giving an improvement of A% 2 « 110. The model parameters are 
quoted in Table [3j The four models are plotted in Figure H] upon the background subtracted 
count rate spectrum of MR 2251—178 in the wavelength space (re-binned to have at least 
100 counts from the source per bin). 



3.3. A Dusty Warm Absorber? 



Dusty warm abs orbers have been found t o imprint their hallmar k on the X-ray sp ectra 



of several AGN (e.g. iKomossa fc Finkl 1 199 71 : Komossa & Bade! Il998l : iLee et al.l l200lh . We 



investigate the possibility of the presence of a dusty warm absorber by looking at the position 
of two important edges, predicted to be strong and noticeable in the X-ray band; the neutral 
carbon edge at 291 eV (from the graphite species), and the O I edge at 538 eV (from silicates). 
Figure [5] shows the position of these two edges (the edges corresponding to the H- and He-like 
oxygen features are included in the model but they lie off the figure axis), in the observer's 
frame (also marked is the Galactic neutral oxygen as O I ga i). These edges are not clearly 
seen in the spectrum, but this does not preclude their potential addition and improvement 
to the model in the fitting. In Table H] we present the result of fitting the spectrum of 
MR 2251—178 with a model that is composed of a Galactic absorbed power-law modified by 
four edges, the K-edges of C I, O I, O VII and O vm. In order to search for fine features, 
we used a larger number of bins in the fit, 50 counts from the source per bin. In terms of 
goodness of the fit, this is equal to the pi x 2 edges model, with xl = 1-22 (for 587 d.o.f). 
However, letting the four edges vary freely, the C I optical depths has no lower limits. An 
upper limit can be set on the dusty material if it exists. We have enforced the tq \ to be 
equal to 0.08 (the upper limit of this measurement), letting all the other optical depths free 
to vary, and re-fit. We set this as an upper limit on the tq \ parameter, since the fit does not 
get worse with respect to the previous one (see row 2 in Table H]). This translates in a C I 
column density Ac i ~ 2 x 10 19 cm~ 2 . Doing the same with the To i (taking the upper limit 
of the best-fit), we obtain no significant change with respect to the best-fit x 2 ( see row 3), 
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and we consider tq i = 0.09 as the upper limit on tq i, which translates in Ng [»9x 10 



19 



cm -2 (we used the photoionization cross sections from iMorrison &: McCammonl Il983l . at 
their respective energies). For typical Galactic gas/dust ratios, our observed upper limit is 
not sensitive enough to provide significant constraints on the dust in MR 2251—178. 



3.4. Mean Warm Absorber 



In this section we discuss the global view of the spectrum of MR 2251 — 178, in terms 
of a more physical model. Visually, the strongest spectral feature is the broad absorption 
structure around 1 7 A (see Figure ED, recognized as the hallmark of the warm absorber (e.g., 
George et all 119981 ) . Detailed calculations have demonstrated that this structure is made 
(mainly) by the O yil-K edge and by the Fe M-shell 2p— 3d unresolved transition array (UTA) 
between 16—17 A ( Behar et al.ll200ll ). The ionization condition under this structure is strong 
enough to dominate the e mergent spectrum, have been the sub ject of several theoretical and 
observational works (e.g.. lKrongold et al.ll2003uGu et al.ll2006l ). due to in part, the particular 
sensitivity of this feature to the ionization state of the gas responsible for the absorption. 

To describe the state of the gas, we build a grid of photoionization models assuming 
ionization and thermal equilibrium. Under conditions of ionization equilibrium the state of 
the gas depends mostly (apart from n H , abundances and column density) up on the shape 



47rF ir 



of the ionizing spectrum and the ionization parameter £, that we define as in iTarter et al. 
jl969h : 

c — (1 

where F ion is the total ionizing flux (F ion = L ion /[47rr 2 ], see below for definition of L ion ), and 
n H is thegas density. We have carried out all the photoion i zation calculations using the 
XSTAR I 11 ! code with the atomic data of iBautista &: Kallmanl ( 20011 ) . The code includes all 
the relevant atomic processes (including inner shell processes) and computes the emissivities 
and optical depths of the most prominent X-ray and UV lines identified in AGN spectra. 
Our models are based on spherical shells illuminated by a point-like X-ray continuum source. 
The input parameters are the source spectrum, the gas composition, the gas density %> the 
column density and the ionization parameter. The source spectrum is described by the 



11 Version 2.1kn6. UTA s included . Howe ver, this version does not include the corrections (for some atomic 
line transitions) given bv lGu et alj (|20061 ). We have re-computed our models bas ed on an updat ed version 
of XSTAR (kindly provided by T. Kallman) which does include the corrections of iGu et alJ ( 20061 ) . We find 
that the change in the outflow velocity is only ~ 1 %. Part of the rea son is that the LETGS has a spectral 
resolution of « 50 mA, while the changes reported by Gu et al. ( 20061 ) were ~ 15 — 45 mA. 
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spectral luminosity L e = L ion / e , where L ion is the integrated luminosity from 1 to 10 00 Ryd, 



and / 



1000 Ryd 



f e de = 1. The spectral function is taken to be the ionizing spectrum of iLeighly 



(120041 ) . The gas consists of the fo llowing element s H, H e, C, N, O, Ne, Mg, Si, S, Ar, Ca 
and Fe. We use the abundances of iGrevesse et al.l (119961 ) in all our models (we use the term 
solar for these abundances). We adopt a turbulent velocity of 100 km s -1 , and a gas density 
% = 10 8 cm~ 3 . 



3.4-1- Global Absorption 

We have fit the X-ray spectrum of MR 2251 — 178 with our XSTAR based photoionization 
models. First of all, we find that independent of column density, a relatively low ionization 
parameter is needed; else, the spectrum in the Fe UTA band is not well reproduced. We 
demonstrate this point in Figure [61 It shows the resulting absorption spectrum for several 
ionization states. For illustrative purposes, here we fix Ah = 4 x 10 21 cm -2 . In the top panel 
we can see the spectrum resulting at high ionization parameter, log(£) = 1.5 and 2. The 
bottom panel shows the spectrum at lower ionization states with log(£) = 0.5, 0.7 and 0.9, in 
black, red and dark blue respec tively. From this w e can rule out high ionization states alone 



to describe this band (though iKaspi et al.l |2004j, reported the use of a multi-component 



model to represent the spectrum of MR 2251 — 178, we were not able to find a significant 
improvement on the fit over the one-absorber component), because these models always give 
large residuals in the UTA band (see below). 

Nevertheless, for our best-fit, single component photoionization model (pi x xstar model = 
Xi, Figure [7]), we obtain the following parameters: a column density of Nn = 4.771^29 x 10 21 
cm -2 , an ionization parameter log(£) = 1.83loo6 an d a photon index T x = 1.96+qq 2 . While 
the fit well describes the long wavelength part of the absorption feature (with respect to the 
centroid around ~ 17.4 A), there exists a large residual in the shorter wavelength range (> 
50 %), resulting in a fit with xl = 1-54 (for 587 d.o.f). Motivated by the good result of 
including a thermal component to describe any soft excess, and also by the fact that the 
X-ray spect rum of this object has been fitted before including a thermal component (for 



example in iGibson et al.l 120051 ) . we include in our model a blackbody, characterized by its 
temperature kT (i.e., (pl+bb) x xstar mo d e i = X 2 ). The inclusion of the blackbody im- 
proved the fit, which gives xl — 1-27 (and d.o.f = 585). But we went further, and inspected 
the possibility that the absorbing gas in MR 2251—178 is flowing outward from the center 
of the system with velocities ~ — 2000 km s _1 . For that purpose we shift the modeled 
spectrum (the absorber) by a grid of velocities keeping fixed the outflow velocity and letting 
the others parameters of interest free to vary. For Xi we found that it favors an outflow 
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velocity 200 km s _1 . In the case of X2, the model points to an outflow velocity ~ 1200 
km s _1 . After finding statistically motivated outflow velocities for our models we were able 
of thawing them and find the best-fit outflow velocity. 

The xstar mo dei plus the blackbody (X 4 ) gives a xl — 1-22 (with d.o.f = 584), and the 
following best-fit parameters: N H = 2.25±g;|| x 10 21 cm" 2 , log(0 = 0.631^6, kT = 89.5±f;| 
eV, and v out = — 1100*2io km s_1 - An inspection of Figure [7] reveals a good agreement 
between the data and the model, and the removal of the residual around ~ 17 A (see ratio 
of data/model) 0. Table summarizes the parameter values of the model with and without 
the thermal component, at rest (i.e., km s _1 ) and with the best-fit outflow velocity. We 
present in Figure [8] a high resolution (bin size 25 ml) version of our best-photoionization 
model plotted over the spectrum of MR 2251 — 178 (in the most interesting range to search for 
spectral lines, ~ 1 —40 A). Here we fix the velocity outflow to the best-fit value found above 
(i.e., —1100 km s _1 ) and allowed the others parameters free to vary (i.e., T x = 1.64^ ; 2, 
N H = 1.82jft| x 10 21 cm" 2 , log(f) = 0.57t° % kT = 95.2±|| eV). We use our best-fit 
featureless continuum (but including bound-free transitions) as the underlying continuum 
for the search for possible spectral lines (also plotted in Figure EJ), discussed in the next 
section. 



3.4-2. Absorption lines 

There are several coincidences between the theoretical prediction about where a line is 
located, and deviations of the data from the fit continuum (the XSTAR featureless global 
continuum described above). In order to seek for narrow absorption features (in a systematic 
way) in the spectrum, we use a finer binned spectrum (bin size 25 mA). We use the following; 
detection criteria: 

1. The counts present in the data must deviate from the fit continuum by at least la. 
As underlying continuum, we use the best-fit global featureless continuum, described 
above. We use a global continuum fit because in that way we are less sensitive to local 
noise in the spectrum 0. 



12 In order to see if a multicomponent warm absorber system is supported by the data we repeat the above 
procedure adding a second warm absorber system to the models, without (X5) and with the blackbody (Xg). 
In general the data do not well constrain the parameters of the second component. And in any case, no 
significant improvement is seen by adding a second component to our model (i.e, \ 2 ~ 707, A\ 2 ~ 6 better 
than X4), and we do not discuss this model further. 

13 As a test of the robustness of our results, we have also determined the continuum locally around the 
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2. The FWHM (within errors) of the line should not be significantly smaller than the 
instrumental FWHM, at the correspondent wavelength. 

Not as a criterion, but as an extra piece of statistical evidence for the presence of lines, 
we use the maximum likelihood ratio test (MLR) to compare the continuum model with the 
continuum+line model, and compute the line significance (see Table |6] column 4). P(> T) is 
the probability that one would select the more complex model (continuum+line) when in fact 
the null hypothesis is correct (continuum alone). We set a lower threshold of P(> T) < 0.01 
(not accounting for number of trials). In fac t, the four feature s which passed criteria 1 and 



2 also had P(>T) below this threshold (see lBand et al .1119971 . for justification of using the 
MLR test). 

A list of candidate absorption lines is given in Table El We measure the center (column 
1), and the width (in terms of a [in km s -1 ], column 2) of each spectral structure using the 
best-fit continuum described above. 

Figure [TO] shows different portions of the spectrum containing absorption lines. Among 
these lines, we only report measurements from four lines (Table ED, located in a region where 
the calibration of the LETGS is reliable and its sensitivity is high enough to report 90 % 
confidence limits on their errors. Deviations > la, but have FWHM too low to be reported 
as individual lines are still marked in the figure (cx-panel), but are not discussed here any 
further. 

The measurements of the line parameters along with their identification is given in Table 
El In the first column we have the center of the line, and in the second column the equivalent 
width of the line (EW) in milli-angstroms, computed with: 



EW 



dX, (2) 



where F S (X) = F C (X)(1 — r exp[— %^ ]), F C (X) is the continuum at the wavelength A, r 
is the depth at the center of the line, a is the measured width in units of A, and Ao is the 
wavelength of the line in A (at the core). The EW uncertainties are computed using the lower 
and upper limits of r and a (i.e, EW min = EW[r min , cr min ] and EW max = EW[r max , cr max ]). 
In the fourth column we give the blueshift of the line (in km s" 1 ). All these are resonance 



line. The continuum was fit locally (±0.5 A) as a power-law. We find that all four key line candidates are 
still present when a local continuum is used. To give an example: the equivalent width of the line of 
N vi (see Table EJ) changes from 55^30 mA to w 70 mA, when using a local continuum. The measurements 
of the other three lines are unaffected by the local continuum fit. 
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lines, product of electric dipole transitions with the form Is — np for C and Is 2 — lsnp for 
N ions. 

We discuss the issue of line identification in two step. (1) As a first (simplifying) step 
we assume, that the spectrum is dominated by only one main warm absorber, and that 
the column density and ionization parameter from the global fit to the X-ray spectrum 
characterize this main warm absorber reasonably well. Then, this best-fit global warm 
absorber model guides us in identifying the absorption lines seen in the spectrum. In this 
first step, we only consider line identification with transitions, which are actually being 
predicted to be strong (i.e., detectable) by this photoionization model. If no line feature at 
zero velocity was predicted by the model, we allowed for a range of velocities (up to ~ few 
thousand km/s - a range of outflow velocities commonly observed in known warm absorbers). 
That procedure resulted in the line identification reported in Table [7J To give an example 
of a possible alternative identification, we compare the strength of the line N vi at 24.9 A 
(ls-2p, the one identified in Table [7J with a potential alternative identification, the line N VII 
at 24.8 A. The former line is much stronger than the latter: The ratio of optical depth at 
the core is r(N vi)/r(N vil) m 8, for our best-fit global ionization parameter of log(£) = 0.6. 
At this log(£), there is almost 3.5 times more N +5 than N +6 . 

Assuming the one- absorber model is correct, the predicted strongest lines are the four 
line candidates presented in Table [7J However, we find that lines coming from the same 
ion are at different outflow velocity. This is a problem at the moment, and opens up a 
second possibility, which we examined in a second step: (2) the possibility that we see 
several different warm absorber components each with a different outflow velocity, and each 
characterized by a potentially different column density and ionization parameter. The data 
quality, a single grating exposure of only 80 ks, does not allow to fit such a four-component 
absorber with so many free parameters. Therefore in step (2) we made the (simplifying) 
assumption that the four identified lines are all at the same velocity. We therefore assumed 
that all four lines are at zero velocity or at velocities up to ~ 500 km s _1 , in particular. 
Independent of any warm absorber model (i.e., any specific £ and Nh), we searched line lists 
to find any consistent line identifications for the four lines, fixing their wavelengths to the 
quasar rest-frame (i.e., assuming zero velocity), or allowing velocities up to ~ 500 km s _1 . 
In that case, no line counterparts can be identified. 

We therefore continue with solution (1), and discuss some implications of this solution. 
Since most warm absorbers are not at rest, in order to identify lines, we allowed for a range 
of outflow velocities, from — 5000 km s" 1 . We find that not all of the absorption line 
candidates are at the same outflow velocity, but we preliminary identify 3 velocity systems 
at ~ —600 km s _1 , —2000 km s^ 1 , and ~ —3, 000 km s _1 . Further discussion of these three 
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components, and a comparison with UV observations is provided in £j5l 

4. Variability in the band (0.6 — 1) keV 

The region ~ (0.6 — 1) keV in the spectrum of MR 2251 — 178 is complex, due to the 
presence and blending of three features: the Fe UTA, and the K-edges of oxygen O VII and 
O Vlll. As a first step, we investigate temporal changes assuming the spectrum is actually 
dominated by the O VII and O vm edges. Later in §5.21 we carefully discuss the variability of 
the spectrum based on a photoionization model, which includes a proper treatment of these 
features including the UTA self-consistently. First, we split the observation in four (equally 
distributed, ~ 11 ks of GTI) time intervals (1-4) and applied a pi x 2-edges model to each 
set of data. By fitting this model (with the edges energy fixed at the values of O VII and 
O vm K-edges, 739 and 871 eV respectively), we are able to notice changes (in time) around 
the bound-free features, at ~ 0.69 keV for O VII and 0.81 keV for O vm in the observed 
frame which in turn modify the parameter values of the model (see Table [S]). To quantify 
these changes we looked at the parameters of the model at each time interval. Due to the low 
S/N of the data these measurements have large errors. Nevertheless, one real inconsistency 
is seen in the depth of the O Vlll edge (toviii) from time interval (1) to (3): approximately 
a factor of 3.4 in Toviii leads to a change in the count rate level of ~ 20 %, in the vicinity 
of the edge. Almost no change is seen in 7"ovn> being consistently about the mean best-fit 
value found before; tovii ~ 0.6. 

We caution that the edge model is a too simple representation of the spectrum in the 
(0.6 — 1) keV range, where also the UTAs are located. Therefore, we do not draw any physical 
implications from it. We discuss a more physical photoionization model in §5.21 

5. Discussion 

In this section we discuss the physical implications of our measurements, which allow 
us to draw conclusions about the kinematics of the outflowing gas, ionization state of the 
warm absorber, and variability of the absorption structure around the oxygen edges seen in 
the spectrum of MR 2251-178. 
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5.1. Identification of absorption lines 



Identification of absorption lines in the X-ray spectra of AGN is important because these 
lines carry a wealth of information on the physical conditions in the ionized gas, and its link 
to other components of the active nucleus. 

We have systematically searched for absorption and emission line features in the X-ray 
spectrum of MR 2251—178, and have identified four absorption line candidates correspond- 
ing to transitions in the ions C VI and N vi. These lines, if real, would imply a rather 
complex velocity field of the absorber, with lines from the same ion indicating different out- 
flow velocities. In particular, it is of our concern the lack of corresponding absorption from 
the element OxygerJ^I. In case of solar abundances of the ionized gas, Oxygen absorption 
features are expected and have been detected in several AGN wit h warm absorbers (e.g 



Kaastra et a 



Fields et al 



200fl 
20071 : 



Sako et al 



2001 



Costantini et al. 



Lee et al 



2007 



2001 



Blustin et ah 2003; Nctzer et al. 2003 



Krongold et al.1 120071 ) even though few of them 



with well-measured column densities. On the other hand, it has also been found that warm 
absorbers are complex, stratified, multi-component, time-variable, and possibly of non-solar 
abundances, and the same may well hold for MR 2251—178. With these cautious comments 
in mind, we briefly discuss implications of the line candidates we find. Ultimately, deeper 
grating observations are needed to study the absorption line features and their temporal 
evolution. 

Independent information on the presence and velocity fields of ionized absorption com- 
ponents comes from a comparison of X-ray with with UV observations. Firstly, in our LETG 
spectrum we identified a candidate gas component outflowing at a velocity v ~ —600 km 
s _1 , as indicated by the C VI A33.736 line. This is consistent with UV measurements. The 
UV spectrum of MR 2251—1 78 provides evide nce for an absorbing system at an outflow 
velocity of v — —580 km s _1 (IKaspi et al.ll2004j ). This component is seen through the O VI 
A1032,1038 and Ly Q lines, observed with the STIS on board HST. Also (and reported by 
the same authors), this component is in accordance with the velocity of ~ —600 km s _1 , 
reflected by the Lyg line, observed with FUSE. This is interesting, because i t confirms the 



scena rio where UV and X-ray absorber systems could be physically related (IMathur et al. 



19981 ). However, the centroid of our X-ray line candidate is measured with an uncertainty 
that also allows to relate th is X-ray absorption system with the average outflow velocity of 
300 km s" 1 reported bv lMonier et al.1 feoOlh using the Ly a , N v A1240, and C iv A1549 



14 we have care fully searched at the locati on of O VII a t A21 .8 A, and cannot identify any feature. This 
is consistent with Gibson et al. ( 2005h while Kaspi et al. ( 2004 ) report some Oxygen absorption features at 
- 16 - 22 A in the RGS spectrum of MR 2251-178 (their Table 3), albeit with low significance. 
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lines in the ultraviolet (UV) band (with the HST-Faint Object Spectrograph observation of 
MR 2251-178). 

For the computation of the mass-outflow rate, we assume that the outflowing material 
in MR 2251 — 178 forms a spherical shell (n ~ r 2 ) expanding at velocity v. Using the 
approximation N ~ nr, we can write: 

M out = 4iTH H Nrvf cov , (3) 



or 

" ''out 



0M8N ^Aim^)^ M -^ (4) 

where \iu is the mean mass per H atom (equal to ^j), N 2 i is the column density of the 
material in units of 10 21 cm -2 , rig is the distance of the gas in units of 10 18 cm, and / cov is 
the covering factor (fi/[47r]). Using our best-fit Nn ~ 10 21 * 35 , an outflow velocity v = 1100 
km s -1 , at a distance of 1.5 x 10 18 cm (distance chosen for illustrative purposes) we have a 
mass-loss rate M out ~ 0.18/ cov M Q yr -1 . 

In Table [9] is quoted our estimate of the mass-loss rate of the LETG X-ray warm absorber 



system compared with previous values. To compute the mass-loss rate, iGibson et all (120051 ) 
used a high- velocity outflow v ~ 13,000 km s _1 , concluding that M out could significantly 
exceed the accretion rate (~ 0.2 M yr -1 ), by ~ 1 — 2 orders of magnitude, even tak ing the 



Fe xxvi line (see Table 7 of that paper). On the other hand, iMonier et all (120011 ). using 
an outflow velocity of 300 km s™ 1 , estimate a mass-loss rate of ~ 0.9 M yr -1 for a 10% 
covering factor, putting this value close to the accretion rate. Our computation of M® is 
more in agreement with the value found by these last authors. The kinematic energy rate 
carried away by the flow spans ~ (0.07 — 1) X 10 42 / COV erg s -1 . 



5.2. Absorption variability during the observation 

Warm absorbers are known to vary on short and long time scales. We find indica- 
tions for variability during the LETG observation of MR 2251—178 in the wavelength range 
that includes the oxy gen absorption edges and the Fe UTA features. Both, the UTAs (e.g., 



Krongold et al.ll2005l ) and/or the absorption edges could be variable. Mechanisms of vari- 
ability include changes of the ionization state of the absorber to changes in the ionizing 
continuum, changes of the internal structure of the absorber, and clouds crossing our line of 
sight. 

In order to find clues on the variability mechanism, we have fit our most successful warm 
absorber photoionization model (model X 4 ; i.e., (pl+bb) x xstar, with u out = —1100 km 
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s -1 ) separately to the four subsets of the total observation (Fig. [2]). We have explored 
the effects of variable column density, ionization parameter, and intensity of the continuum 
(black body and powerlaw) spectral components. The strongest effect is a change in column 
density of the ionizing material in epoch (3); formally requiring (at > 6a) a smaller column 
density in order to fit successfully that epoch (Fig. [TTI) . In order to re-check whether the 
data do indicate variability between the different epochs, we have used the best-fit model of 
epoch (1) and compared it to the data of epoch (3). In the comparison, we have fixed all the 
parameters (of our model X 4 ) to those derived for epoch (1), and then inspected the residuals 
during epoch (3). Clear deviations are seen (Figure [T2l). demonstrating independently the 
presence of spectral changes throughout the observation. 

Changes in the luminosity would affect the ionization state of the absorber, but the 
average countrate or MR2251-178 does not vary throughout the observation. Taken at face 
value, a true change in column density would require a change in density and/or thickness 
of the absorber, and we first briefly comment on this possibility. Assuming the thickness 
of the absorber has not changed in such a short time, the density of the absorber would 
need to change by a factor of 3 to account for the change in Nh between periods (1) and 
(3). However, changing the density of the whole extended gas (or a fraction of the gas by 
a large amount), is not easy in such a short time scale. Re-arrangements of clouds crossing 
our line-of-sight would have to occur within ~ 10 ks, which is very unlikely. Alternatively 
we may see non-equilibrium ionization effects in the gas, in which the material responds 
with some time delay in a complicated way to previous variability in the ionizing continuum, 
which was not within our observation time window. Deeper observations with longer time 
bases are required to investigate these possibilities further. 



6. Conclusions 

The measured (0.1 — 2) keV flux /(0.1-2 keV) = 2.58+o;q 4 x 10~ n ergs cm -2 s" 1 implies 
that MR 2251—178 is in a relatively low state. The soft X-ray luminosity amounts to 
£(0.1-2 kcv) = 2.54 x 10 44 erg s -1 . 

We did not find any strong evidence for the presence of an extra cold material with 
column density of Nh ~ 10 21 cm -2 toward MR 2251—178 (if part of this is located along our 
line of sight). Based on different spectral fits, we set the upper limit of this component to 
iV co id « 1.2 x 10 20 cm" 2 . 

As for previous observations of MR 2251 — 178, power-law plus blackbody does not pro- 
vide a successful fit to the X-ray spectrum. The addition of a warm absorber improves the 
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fit significatively. 

Based on XSTAR photoionization modelling, we find the observation to be consistent 
with ionized absorbing material in our line of sight, column density iV H ~ 2 x 10 21 cm -2 , 
and an ionization parameter of log(£) ~ 0.6. The inclusion of an additional thermal compo- 
nent is important because otherwise: (1) the absorption structure around 18 A is not well 
reproduced, (2) the ionization state of the gas would be higher, bringing a lot of problems 
in the identification of the atomic line transition. The temperature of the blackbody used 
to represent this thermal component is kT ~ 90 eV. 

We find four line candidates at moderate confidence level. If these lines are real, their 
presence would imply three components traveling at velocities ~ 600, 2000 and 3000 km s~ 4 . 
We compute M ~ 0.01 — O.1M yr _1 and a kinematic energy of ~ 1 o 40-42 erg s _1 , using the 
C VI and N VI lines outflowing at velocities of ~ 2000 — 3000 km s -1 . Using f cov = 0.1, this 
value is a factor of ~ 20 less than the accretion rate of the system. However, we caution that 
these candidates need to be confirmed with future observations. The lack of corresponding 
oxygen absorption is a problem at present. 

We do not find positive evidence for a dusty warm absorber. The spectrum of MR 2251 — 178 
allows us to set an upper limit of N c r « 2 x 10 19 cm -2 , and N Q i ~ 9 x 10 19 cm -2 , if dusty 
material is present in the nucleus of MR 2251 — 178. 

We find changes in the absorption structure of MR 2251 — 178 during the observation, 
not accompanied by changes in the observed continuum luminosity. Possibly, we see non- 
equilibrium effects in the ionization of the gas responding to previous changes in luminosity, 
or changes in density; the true mechanism can only be uncovered with deeper follow-up 
observations. 

We thank the anonymous referee for many constructive comments. The observation 
would not have been possible to analyze without the enormous effort of the Chandra team. 
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Table 1. Chandra Observation Log of MR 2251-178. 



Sequence Number UT Start UT End Time (ks) 

700405 2002 Dec 23, 16:35 2002 Dec 24, 14:54 78.5 



Table 2. Cold absorption fit results 



A H a T x Norm 11 T x (2) Norm h (2) kT 1 Norm s x£/(d.o.f) 



1 fi8+°- 07 b 
l.D8_ 007 


1 70+0-02 
L ' a -0.02 


r o 4 +0.06 








2.55/306 


2.77 c 


i 04+0-02 

i - y ^-0.02 


J - UJ -0.06 








3.00/307 


3.27 d 


2 oo +0 - 02 

z - uu -0.02 


t- 70+O.O6 
°- '°-0.06 








3.40/307 


3.77 c 


z - ud -0.02 


5 qi+0.06 

°- y± -0.06 








3.85/307 


< io- 6 b 


i p-9+0.02 
1 - dz '-0.02 


4 7n+0-0 7 

+• ' u -0.05 






••• 80.31?;° 1.0-1!!- 


1.64/304 


i 9/1 +0.10 b 

l - zq: -om 


o 47+0.03 
-0.03 


i 20 +0 - 03 

1 -^ u -0.03 


1 


a-i +0.03 
4i -0.03 


^ 07+0-06 
°- y ' -0.06 


1.81/304 


1.74 f 


o 7n+0.03 
°- ' u -0.03 




1 


44+0-03 

^-0.02 


4 1 fi+006 
^• iu -0.06 


1.81/305 


2.24 g 


o n 4 +0.03 
°- y ^-0.03 


o qi+ - 02 

u - y± -0.02 


1 


47+0.03 
^' -0.02 


a q/i+0.06 
0.06 


1.84/305 



Note. The error parameters are 90 % confidence limits, (a) Absorber column 

density in units of 10 20 cm -2 , (b) Free parameter, (c) Fixed to Galactic value, (d) Fixed 
to Galactic value plus A Ah, where AAh = 5 x 10 19 cm -2 , (e) Fixed to Galactic value plus 
2 xAA+i. (f) Extra absorber column density (shifted by z = 0.06398, model zphabs in 
XSPEC), fixed to the best value plus AA H . The fit is corrected by the Galactic absorption, 
(g) Extra absorber column density (shifted by z = 0.06398, model zphabs in XSPEC), 
fixed to the best value plus 2 xAA H . The fit is corrected by the Galactic absorption, (h) 
Powerlaw normalization, xl0~ 3 photons keV _1 cm" 2 s _1 at 1 keV. (i) Temperature of the 
blackbody in eV. (j) xl0~ 6 in units of L 39 /Dl , where L 3g is the luminosity of the source 
in units of 10 39 erg s" 1 , and D w the distance to the source in units of 10 kpc. 



Table 3. Results from simple model fits 



Model b 


r x 


Norm c 


r O VII 


TO VIII 


kT d 


Norm 6 


Xl /(d.oi) 


[1] 
[2] 
[3] 
[4] 


1 04+0-02 
i - y ^-0.02 

1 [-0+0.02 

i - d °-0.02 

1 Q9+ - 02 
i - oo -0.02 


c «^+0.06 
d -°°-0.06 
70+0. 06 
^•'°-0.06 

u -°°-0.07 
fi ^ +a07 


u -°°-0.04 
u - di -0.04 


n 9 9 +0.05 
u - zz -0.04 

21 +0 - 05 


81.5115 

44.1118 


y - y ^-0.38 

cr p-7+0.92 
°-° 1 -0.92 


3.00/(307) 
1.63/(305) 
1.38/(305) 
1.25/(303) 



Note. All the models include absorption fixed at the Galactic value of Nn = 
2.77 x 10 20 cm -2 . The error parameters are 90 % confidence limits, (b): Model [1]: pi; 
[2]: pl+bb; [3]: plx2 edges; [4]: (pl+bb)x2 edges; pl=powerlaw; bb=blackbody. 
(c) Powerlaw normalization, in units of 10~ 3 photons keV _1 cm _2 s _1 at 1 keV. (d) 
Temperature of the blackbody in eV. (e) xlO -5 in units of L 39 /Df , where L 39 is the 
luminosity of the source in units of 10 39 erg s _1 , and D w the distance to the source in 
units of 10 kpc. 
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Table 4. Dusty warm absorber fit results 



r x 


Norm b 


T"0 VII 


TO VIII 


TC I 


TO I 


xl/(d.osy 


1 Q9+ - 02 

1 QQ+0.02 
± - yo -0.02 
1 OQ+0.02 
1 - yo -0.02 


v -°°-0.07 
u - do -0.07 

fi ss +0 07 


u - od -0.04 

o.55ig;8i 
o.54i8:8^ 


n oo+O. 04 
u - zo -0.05 
n oo+O. 05 
u - zo -0.04 

24 +0 - 05 

u - z ^-0.04 


0.02+o.oe 

0.06 c 
004+ o.o6 


0.0518:81 
0.0618:81 

0.09 d 


1.22/585 
1.22/586 
1.22/586 



Note. - The error parameters are 90 % confidence limits, (b) Powerlaw 
normalization, xl0~ 3 photons keV _1 cm _2 s _1 at 1 keV. (c) Fitting, enforcing 
r c i = 0.08, and r i free, (d) Fitting, enforcing To i = 0.09, and r c i free, (e) 
In order to search for fine features, we used a larger number of bins in the fit, 
50 counts from the source per bin. For comparison the best fit model of the 
pi x 2edges model (using this data) gives xl — l-22/(d.o.f=587). 
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Table 5. Fit results for photoionization models. 



Parameter Xi X2 X3 X4 



r x 


l-96i°;°J 


1 66+ 002 


1 07+0.02 
1 - y ' -0.02 


1 66+ - 02 

i - uu -0.02 


Norm a 


6.84l8g 


r 42 +0.07 


6.821H? 


c qo+0.07 
d - oo -0.07 


kT (eV) 




89-9±i:I 




89.5±1;1 


Norm bb b 




in 21 +0 - 42 




in.29l^ 


iV H c 


4 77+0.31 

' ' -0.29 


2 44+ - 13 

z -^-0.29 


4.28_° ) ; 2 ? 


9 9 t-+0.28 
z - zd -0.25 


log(0 


1 00+O.O6 
± - oo -0.06 


n 64 +0 - n 


1 7 9+0.07 
' z -0.07 


0.63_8£ 


^out d 





n 


-34 _ 2 20 


-llOOlgo 


^/(d-o.f) 


1.54/(587) 


1.27/(585) 


1.52/(586) 


1.22/(584) 



Note. - (a) Powerlaw normalization, in units of xin~ 3 pho- 
tons keV~ 1 cm" 2 s _1 at 1 keV. (b) xin~ 5 in units of L 39 /Dj , 
where L 39 is the luminosity of the source in units of in 39 erg s^ 1 , 
and D w the distance to the source in units of in kpc. (c) Col- 
umn density of the ionized material in units of xin 21 cm~ 2 . (d) 
Outflow velocity of the ionized gas in units of km s _1 . Models 
Xi — X 4 are described in the text (see section [3.4. ip . 
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Table 6. Gaussian absorption line candidates 



Obs. A 

(A) 


^measured 

(km s _1 ) 


C"LETG 

(km s _1 ) 


MLR 

P(> T) 




26.341°;°? 
30.091°;°* 

on qo+0.03 

OU - JZ -0 .02 

ok OO+0.03 
OO.OO_ 03 


2801?" 
22011?° 


242 
212 
210 
178 


7.38 x 10" 
1.13 x 10~ 
7.06 x 10" 
1.17 x 10" 


-3 
-3 
-4 
-2 



Note. Error parameters are 90 % confi- 
dence limits. The final significance is less than 
the one shown in column 4, since the number 
of trials is not considered. The MLR test is de- 
scribed in Section 13.4.21 
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Table 7. Identification of absorption line candidates 



^obs 


EW 


Ion and atomic 


Alab 


Vshift 


(A) 


(mA) 


transition 


(A) 


(km s _1 ) 



26.3418;" 


°°-30 


N vi Is 2 X S - 


IsSp 1 ? 


24.914 


-20301^8 


30.091^2 


63^ 


C vi Is 2 S 


- 3p 2 P° 


28.466 


-2080l 2 }° 


on qo+0.03 

ou.oz_ 02 




N vi Is 2 X S - 


1S2P 1 ? 


28.787 


-3220 + ^° 


35.83i° : °i 


86l|g 


C vi Is 2 S 


- 2p 2 P° 


33.736 


-580±I?g 



Note. - 



Errors are 90% confidence limits. 



Table 8. Variability of the oxygen edges" 



Period 


''start 


''end 


total 


r x 


Norm c 


r +e 


Tq+7 


xl/(d.ol) 




xlO 8 (s) 


xlO 8 (s) 


counts (xlO 4 ) 













(1) 1.57048747 1.570653992 1.8298 

(2) 1.570653993 1.570819941 1.6104 

(3) 1.570819942 1.570982180 1.7391 

(4) 1.570982181 1.571276092 1.5842 



1.90l°°4 6.11+°^ 0.57t°"°9 0.09+ - 08 1.19/161 

i.89igSJ 6.58+8:i! o.67±8:5S 0.25^ 1.01/145 

1.9118$ 6.621°'^ 0.461°-™ 0.3ll°™ 1.32/159 

i.96i^ 6.6ii8:2 0.5718:^8 0.221™ 1.39/144 



a The model includes Galactic absorption of Nn = 2.77 x 10 20 cm 2 . The error parameters are 90 % confidence 
limits computed with the spectrum adaptively binned to have at least 50 counts from the source per bin. 

b In the detector (original) time. These intervals exclude data from the "bad time intervals" (section 12.11) 
c Powerlaw normalization, in units of 10~ 3 photons keV _1 cm~ 2 s~ 1 at 1 keV. 
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Table 9. Outflow masses and kinematic energy rates 



Ion a 


MMonier et al. 2001) 6 




\Mv 2e 


M f 


\Mv 2i 




[M Q yr" 1 ] 


[M e yr- 1 ] 


[erg s" 1 ] 


[Me yr- 1 ] 


[erg s -1 ] 


C IV 


0.9 










Fe xxvi 




220/ COV 


1.1 X 10 46 /cov 






Fe xvn 




8900/ COV 


7.9 x 10 47 / COV 






S xiv 




1700/ cov 


1.6 x 10 47 / COV 






(c) 








0.18/ cov 


6.7 x 10 40 /cov 


N vi d 








0.32/ COV 


4.0 x 10 41 / cov 


C vi d 








0.48/ COV 


1.4 X 10 42 /cov 



Note. - - (a) Ion from which the rates are computed, (b) For a covering factor of 10%. (c) 
Best-fit outflow velocity found in §3.41 of 1100 km s _1 . (d) Using the velocities shown by the 



10 



21.35 



cm' 2 ) , of ps 2000 km s -1 , and 3000 km s _1 for N, and C 



lines in Table [7] (and iV H 
line respectively, (e) From I Gibson et al.l (120051 ) . (f) Present work. 
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Fig. 1— 0.1-10 keV LETGS spectrum of MR 2251-178 (observed frame). It is binned to 
have at least 100 counts from the source per bin. In red the model of a simple absorbed 
power-law modified by 2 oxygen absorption edges. 
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Fig. 2. — Light curves of the source and background (bin size 200 s). The top panel shows, 
in red, the mean counts of the source. No significant time variability is seen within ±ler 
(green lines). The bottom panel shows strong variability of the background (here, the green 
lines show ±3<r). The good time intervals are indicated by the interrupted line below the 
background spectrum. The meaning of intervals (l)-(4) is explained in §H 
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0.1 1 10 

Observed Energy (keV) 

Fig. 3.— LETGS spectrum of MR 2251—178 (positive and negative orders added) in 
energy space. Top panel: The solid red line shows the absorbed power-law fit. The solid 
blue line is the absorbed pl+bb model. Middle panel: The solid line is the absorbed two 
power-law model. Bottom panels: Ratios of the composite models to the data: absorbed 
power-law (black) and blackbody plus power-law (blue), in the bottom panel the ratio for 
the two power-laws model (absorbed, means cold absorption with zphabs ). 
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Fig. 4. — LETGS spectrum of MR 2251 — 178 in the observed frame (positive and negative 
orders added). It is binned to have at least 100 counts from the source per bin. All the 
models include the Galactic absorption toward MR 2251—178 (see text). Top panel: Solid 
lines: pi (red), pi x 2 edges (green), pl+bb (dark blue), (pl+bb) x 2 edges (light blue). 
The rest of the panels are the ratios of each composite model to the data, represented by 
their corresponding colors (pl=power-law; bb=blackbody). 
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Fig. 5. — Test for the presence of a dusty warm absorber. X-ray spectrum of MR 2251—178 
(observed frame), with a power- law model modified by four edges (solid red line). The edges 
corresponding to the H- and He-like oxygen features are included in the model but they lie 
off the figure axis. Marked are edges expected from dusty material represented by bound-free 
transitions of neutral carbon and oxygen, in the observer's frame. Also marked is the oxygen 
edge coming from the Galactic absorption toward MR 2251 — 178. 
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C3 




10 20 30 40 

MA) 

Fig. 6.— 10 — 40 A spectrum of MR 2251 — 178 and XSTAR photoionization models of 
different ionization parameters. We use Nh = 4 x 10 21 cm~ 2 for illustrative purposes. Top 
panel: High ionization state, log(£) = 1.5,2; black and red lines respectively. Bottom 
panel: Lower ionization state, log(£) = 0.5,0.7,0.9; black, red and blue lines respectively. 
It can be observed that at higher log(£), the produced absorption feature is at bluer wave- 
lengths. 
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Fig. 7. — The best-fit XSTAR photoionization model (blue). A residual is seen at ~ 17 A. 
An important improvement is achieved by the inclusion of a thermal component (red). See 
Table [5] for parameter values. 
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Fig. 8. — Comparison of the spectrum of MR 2251 — 178 with our best-fit XSTAR photoion- 
ization model which is presented in high resolution (bin size 25 mA). The physical model is 
shifted by —1100 km s _1 . Predicted lines of Ne, Fe, O, C and N are marked as labels in the 
top of the spectrum. 
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Fig. 9. — Comparison of the spectrum of MR 2251 — 178 with our best-fit XSTAR photoion- 
ization model which is presented in high resolution (bin size 25 mA). The physical model is 
shifted by —1100 km s _1 . Predicted lines of Ne, Fe, O, C and N are marked as labels in the 
top of the spectrum. 
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Fig. 10. — Strongest candidate absorption lines. The K a and Kp lines of C VI and N VI 
show at least three outflow components at ~ —600, —2000 and —3000 km s -1 . In red the 
Gaussian fit, in blue the XSTAR model shifted by —1100 km s _1 . 
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Fig. 11. — Time variability of the parameters for our model X4 with the column density Nu 
of the warm absorber free to vary, (a) Column density in units of 10 21 cm -2 , (b) Ionization 
parameter log(£), with £ in units of erg cm s^ 1 . 
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Fig. 12. — Best-fit model (X4) resulting from epoch (1) with xt — 11 (interrupted-line), 
applied to the data of epoch (3) (solid-line). Deviations are seen in the soft X-ray band 
(Xl = 1.7). 



